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The string inspired higher dimensional theories suggest modification of Newton's law at submil- 
limeter length scales. Inter-particle distances in white dwarf and neutron stars are 10 -10 cms and 
10 -13 cms respectively, and therefore, the effects of of short distance corrections to gravity deserve 
investigation. We show, by carrying out explicit analytical many-body calculations that, in the pres- 
ence of corrections, the normal state of these compact stars become metastable. The actual quantum 
mechanical ground state of the stars turns out to be unstable. However, the tunneling probability 
to the unstable ground state is so small that the stars may remain trapped in the metastable state 
practically for ever. 

PACS numbers: 04.50.+h 
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I. INTRODUCTION 

The brane-world scenario proposed by Randall and 
Sundrum 0, Q, gives rise to modification of gravity 
at submillimeter distance scales; the corrections are in 
higher powers of inverse distance and appear as addi- 
tional terms in the Newtonian potential. The modified 
potential is given by, 
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(1) 



where I is the curvature scale of 5-dimensional anti- 
deSitter space time, and as before, G is the usual New- 
ton's constant of gravity, M is the mass, r is the distance 
in 3-space. In the setting above, I has the dimension of 
length in 3-space. Henceforth, we will take l 2 s = |Z 2 . 

The astrophysical consequences of short distance mod- 
ifications of gravity are discussed in Refs.0, The 
average inter-particle distance in compact stars such as 



the white dwarf and neutron stars are 10" 



-10 



and 



10 -13 cms respectively, and therefore, the submillimeter 
corrections, in principle, would effect such parameters as 
mass and radius, and also the mass-radius relationship 
of the stars. Wiseman |4| has carried out a detailed nu- 
merical study of the neutron stars on the brane and has 
claimed that these parameters of the neutron stars are 
not effected by brane world corrections. 

In Ref.pl, we have carried out some simple analytical 
many-body calculations using the brane modified grav- 
itational potential. Our approach was based on many- 
body theory in semi-relativistic quantum mechanics. It 
was demonstrated that the inclusion of the brane cor- 
rection terms does not effect the Chandrasekhar ground 
state which is consistent with the results obtained by 
Wiseman using numerical technique In this paper, 
we intend to show that, although correction terms do 
not effect the variational minimum of the energy corre- 
sponding to Chandrasekhar ground state, they give rise 



to two new features of energy as a function of the radius: 
a positive maximum very close to the origin and a new 
unbounded minimum at the origin itself. The unbounded 
minimum indicates the absence of true quantum mechan- 
ical ground state and, there by, a quantum mechanical 
instability. This implies that the normal Chandrasekhar 
ground state of these compact stars become metastable. 
We show that the tunneling probability of the stars from 
the metastable state to the unstable ground state is van- 
ishingly small. 



II. THE KINETIC ENERGY OF RELATIVISTIC 
DEGENERATE ELECTRON GAS 

We consider quantum mechanics of a self-gravitating 
system of large number of fermions in a compact spher- 
ical volume, V, with a large constant density such that 
fermions are relativistic. We will, first, consider here the 
white dwarf stars and, then, argue that our results can 
also be carried over to neutron stars. Consider a white 
dwarf star whose constituents are atoms of mass number 
A and atomic number Z. Let the number of atoms in the 
star be Na- The star is charge neutral, and therefore, the 
total number of electrons in star, N e = ZNa and the to- 
tal number of nucleons (protons+neutrons) N n = AN a- 
Latter, we will see that, N e /N n — Z/A = Y e is an im- 
portant characteristic of a white dwarf star. Our basic 
assumption are: 

(1) For white dwarf stars general relativistic effects can 
be neglected. 

(1) Atomic nuclei are practically motionless and provide 
the entire gravitational attraction. 

(2) The star is so dense that electrons constitute a de- 
generate relativistic Fermi gas at zero temperature. 
The kinetic energy of the star is the energy of the rela- 
tivistic degenerate free Fermi gas which can be evaluated 
as follows 71. The single-particle electron states are spec- 
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ified by the momentum p and the spin quantum numbers 
s = ±1/2. The single-particle energy levels are indepen- 
dent of s : 



e ps = y/ (pc) 2 + (m e c 2 ) 2 - m e c 2 

where m e is the electron mass. The energy of free Fermi 
electron gas is given by, 



£ F = £ £(x F ) - N e m e c 2 ; £ 
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where xp — pF/m e c and the Fermi momentum is given 
by p F = h(3-K 2 n e y/ 3 (n e = N e /V). 

In terms of xp, £q = Sir 2 N e m e c 2 / xp. Therefore, the 
kinetic energy of the electron gas for constant electron 
density is given by 



37r 2 N e m e c 2 
K.hj. — £p = o £(xf) 
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For variable electron density, a similar expression has 
been derived by Daubechies in a very general and mathe- 
matically rigorous way (§, see also [{J). We will mainly 
be considering spherically symmetric stars with a con- 
stant density profile. Following reference jlOj . we express 
the dimensionless Fermi momentum, xf , in terms of mass 
and radius of the white dwarf stars as follows. Let, Mo 
and i? , of dimensions of mass and length be given by, 



Rr, = 
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For helium white dwarf stars, Y e = 1/2, therefore, Mo = 
2.65OM and R = 8623fcm. 

It will also be convenient to introduce the dimen- 
sionless mass and radius of stars as M = M/M and 
R= R/Rq in terms of which xp — M 1 / 3 / R. 



III. THE INTERACTION POTENTIALS AND 
THE POTENTIAL ENERGY 

We consider a spherically symmetric star of mass 
M and radius R, with a constant mass density profile 
Pa{t) = p = constant for < r < R, and Pa(t) — 
for r > R. Therefore, M = N A m A = ^-pR 3 . In what 
follows we briefly describe the potential energy of the sys- 
tem in presence of short range corrections due to Randall- 
Sundrum (RS) and Yukawa corrections (see Ref.@ for 
details) At any point r inside a solid sphere, the poten- 
tial consists of two parts: one part is due to the inner 
sphere bellow the point r and second part is due to the 



outer spherical shell. Let $W(r, a) be the potential out- 
side a solid sphere of radius a, and ^ 2 \r, 6, R) be the 
potential inside a spherical shell of inner radius b and 
outer radius R. The potential energy is given by, 



U =2 



47r / 9r 2 dr($ (1) (r - e,r) + $ (2) (r,r + e,R) 



(4) 



At the end of the calculations, we try to take the limit 
e — > 0. If the limit does not exist, we take e to be the 
smallest inter-particle distance. 

Using the expressions for the potentials <&^g (r, a) and 

$^g(r, 6, i?) obtained in Ref.0], we obtain the potential 
energy corresponding to RS corrections in Eq. (1), 



3 GM 2 l 2 R 
URS = -2-^ ln ~e 



(5) 



Here e is the average inter-particle distance in the star. 
For a system of Fermions, it will be proportional to 
R/N 1 / 3 . The precise numerical factor which is of order 
one, will not matter for our discussions, and therefore, 

1 /3 

we will simply take e = R/N n ■ The potential energy 
can, therefore, be written as 



U, 



3 GM 2 l 2 
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In terms of the variables Mo and Rq, 
U rs = 

Let us call the terms in the first bracket (3. In terms 
of the Fermi momentum, xf the total potential energy 
corresponding brane modified potential given by Eq. (1) 
is 




U = -N P m e c 2 M 2/3 



xp — N e m e c 2 (3x 3 



(7) 



where the first part in the above equation is contributed 
by the Newtonian potential. For submillimeter correc- 
tions, j3 ~ 10 -19 for helium stars. For other stars it 
varies within a order magnitude of this number. 

Yukawa interaction is a very generic short range inter- 
action (see 01 f° r details) and therefore, the results ob- 
tained for this potential is expected to hold qualitatively 
for most short range potentials. On phenomenological 
grounds we require, that the strength of the potential 
due to the correction term at the millimeter length scale, 
at the best, be comparable to gravity. This means that 
the parameter, a, characterizing the strength of poten- 
tial, be a number of the order of unity. We will keep in 
mind this fact while preserving a in all our calculations 

G3 

For spherically symmetric star of radius R and con- 
stant mass density p, the many-body potential at any 
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point r , < r < R is given by fll| . 



1? 



(8) 



Let us note that the second term in the bracket is very 
small and is negligible in the body of the star but it 
is large close to the boundary and will give rise to sur- 
face tension as in the case of Randall-Sundrum correc- 
tions. We are interested in the balance of pressure and 
gravity (along with various corrections) in the body of 
the star, and therefore, we neglect this term. Thus in 
the potential given in Eq.(8), the only relevant part is: 
$y(r) = —4TrGap/fi 2 . Therefore, the potential energy, 



U Y = 



4irr 2 p<^Y{f)dr = —N e m e c 2 jx F 
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where 7 
10 



19 for helium stars 



For submillimeter corrections 7 ~ 
Note that the form of potential 



energy for both RS and Yukawa corrections are cubic in 
Xp. The dimensionless parameters fi (for RS) and 7 (for 
Yukawa) are of same order of magnitude and extremely 
small. Henceforth, we use fi for both. 
Now, the relativistic quantum mechanical total energy is, 
E(M,x F ) = N e m e c 2 E(M,x F ), where 



E(M,x F ) = 



3ir 2 £(x F ) 



1 - M 2/3 



x F — fix' 
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In Eq.(10), the first two terms arise from the kinetic en- 
ergy given in Eq. (2) , the third term is due to the Newto- 
nian potential energy (see reference H3)> an d the last 
term comes, as shown in Eqs.(7) & (9), from poten- 
tial energy due to the short range corrections. In the 
absence of correction terms, the total energy is E c = 
N e m e c 2 E c (M,x F ), where 



E c (M,x F ) 
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IV. GROUND STATE OF THE COMPACT 
STARS IN THE PRESENCE OF SHORT RANGE 
CORRECTION TO GRAVITY 



Eq.(ll), has a unique minimum. As we change the mass, 
the location of the minimum changes. By definition, x F 
is inversely proportional to the radius, and therefore, the 
variation of the location of the minimum with the mass 
gives the mass-radius relationship of the star. 
We choose M = Mi(~ 0.25), such that the star is Chan- 
drasekhar stable with energy minimum at x F = xi(~ 
0.8). Now consider the energy with correction terms, 
Eq.(7), for M = M±. The correction term is too small, 
and become relevant only for very very large values of 
x F and therefore, the local minimum, which exist for 
small x F in the absence of corrections, will persist in the 
presence of the correction terms. Therefore, as before, 
energy as function of x F will have a local minimum at 
x F = x\. However, for very large value of x F , the en- 
ergy approaches —00. To see this let us note, that for 



Xp » 1. 



3^£(x F ) _ 3 



x F . Therefore, 



M 



Xp — 1 — fix 



(12) 



Since Mi < M c , the expression in the square bracket in 
eq.(12) is positive. It is easy to see that E, given above 



has a maximum for x F = 



~ M i /3 )- The value 

of the function E at the point of the maximum is very 
large. Beyond the maximum, E rolls down to —00. To 
summarize, the energy with correction term has the fol- 
lowing structure: for small x F , there is a minimum, then 
for very large value there is a maximum after which it 
approaches —00. The dimensionless Fermi momentum, 
x F , is inversely proportional to the radius, R. Therefore, 
energy as a function of R has a minimum for some fi- 
nite value of R, say R = R\ corresponding to the Chan- 
drasekhar state. In addition it has a maximum very close 
to R = (say, for R = R 2 ), followed by an unbounded 
minimum at the origin. This makes the Chandrasekhar 
state metastable. In the dimensionless unit that we intro- 
duced, i?i is of order one, where as R2 ~ \ffi ~ 10 -9 5 , 
therefore, R\/R% » 1. 

The potential and force due to the short range correc- 
tion terms are similar to those in nuclear physics ex- 
cept that the strength is much weaker and the range is 
much larger. Therefore, a reference to nuclear physics 
may unveil qualitatively the relevant physical features. 
Consider, the energy-per-particle without the Newtonian 
term, E'(M,x F ) = N e m e c 2 E'(M, x F ), where 



A. The white dwarf stars 

Large quantum systems such as white dwarf stars are 
stable if the ground state energy per particle is bounded 
bellow. For a fixed value of M, we look for the extrema 
of per-particle energy, S = E /(N e m e c 2 ), as a function 
of x F . In the absence of corrections, it is known from 
the works of Chandrasekhar see also 0) that bellow 
some critical mass, M c , the energy functional, S c given in 



S'(M, x F ) 



Zii 2 £{x F ) 



1- fix 



(13) 



This function has a maximum and an unbounded mini- 
mum. There is no other minimum. This feature is remi- 
niscent of heavy nuclei, with nucleons interacting via the 
Yukawa potential. With suitable interpretation of the 
parameters, Eq.(13), can, in fact, be carried over to nu- 
clear physics for quantitative treatment of heavy nuclei. 
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It is well known in nuclear physics that the zero-point 
kinetic energy whether relativistic or non-relativistic can 
not balance the attractive potential energy of the Yukawa 
interaction. The problem is solved by invoking a very 
short range hard core repulsive potential which provides 
a lower bound and removes the unbounded minimum. 
Without this hard core repulsive potential, heavy nuclei 
can not be stable. 

In our case we do not have a repulsive hard core poten- 
tial. But the inclusion of long range Newtonian potential 
in Eq.(13) makes it identical to Eq.(lO), and therefore, 
creates, bellow a critical mass, a local minimum in the 
energy function without effecting the other two features: 
a maximum and an unbounded minimum. The original 
Chandrasekhar ground state shows up as a metastable 
state in the local minimum of the energy functional. 



B. The neutron stars 



The neutron stars can be analyzed exactly in the same 
way as the white dwarf stars except for the fact that for 
neutron stars the general relativistic effects would play 
some role. In what follows, we, first, neglect these effects 
and comment on them at the end of this section. The 
first thing that should be noted is that both the kinetic 
and the potential energy in the neutron stars are due to 
the neutrons alone. Therefore, all the formulae and equa- 
tions which are valid for the white dwarf stars are also 
valid for the neutron with the following replacements: 
(electron density, n e — > neutron density, n n ), (electron 
mass m e — ► neutron mass, m n ), (number of electrons, 
N e — * number of neutrons N n ), and (Y e — ► 1). This 
leads to the change in the numerical value of the follow- 
ing quantity, M -> M (n) = 10.599 M , R -> R { Q n) = 
9.40fcm ,/? -> =~ 1CT 13 , 7 -> =~ 1CT 13 . The 
ground state of the neutron stars can be found just as in 
the case of white dwarfs. We reach the similar conclusion 
that, in the presence of short range corrections, neutron 
stars are metastable. 

The general relativistic effects in the compact neu- 
tron stars , in our opinion, influences the Chandrasekhar 
ground state and the values of parameters such as mass 
and radius that are associated with it. It also changes 
the critical mass for the stability of the star. However, 
it should be clear from the discussions in the earlier sec- 
tions that the short range correction to Newton's law of 
gravity does not effect the critical mass for the stability of 
the stars given by Chandrasekhar limit. It also does not 
effect the mass, radius or mass-relationship when the star 
is in the metastable Chandrasekhar ground state. This 
suggests that our results, in all likely hood, are valid in 
the presence of general relativistic effects. 



C. The probability of tunneling to the unstable 
ground state 

Let us provide an estimate of the probability of tun- 
neling from the metastable state to the unstable ground 
state. It is obvious that the peak, at the maximum of 
the energy functional, that separates the metastable state 
(local minimum) from the unstable ground state (global 
minimum), acts as an energy barrier for the tunneling 
process. Let us recall some of the formulae from quan- 
tum mechanics on barrier penetration. Suppose that po- 
tential between two points, x\ and x 2 increases mono- 
tonically, the point, x 2 , is the peak of the potential and 
the energy of the particle crossing the barrier equals to 
its' value at the lower end. Then, the transmission co- 
efficient or the tunneling probability, T ~ e~ s / h , where 
the action, S = f* 2 p.dx, and S/K >> 1. For the com- 
pact stars, let us represent by the radius of the star 
in the metastable ground state and, by R 2 , the radius 
of the star at the local maximum separating the un- 
stable ground state. As argued in the earlier sections, 
R\ > > R 2 . The typical momenta of individual particle 
in these states are, approximately, given by the Fermi 
momenta, hN 1 ^ 3 /R 1 and TiN 1 / 3 /R 2 . Tunneling of the 
compact stars through the energy barrier at the max- 
imum of the energy functional amounts to the change 
in the radius of spheres, occupied by individual parti- 
cles, from Rx/N 1 / 3 to R 2 /N 1 / 3 . This means that inter- 
particle distances change from Ri/N 1 ^ 3 to R 2 /N 1 / 3 , and 
as a result, momenta of individual particles change from 
HN^/Ri to hN 1 ' 3 /R 2 . We will first calculate proba- 
bility for such a tunneling for a single particle. The ac- 
tion S, for this can be given by S = f^ 2 p.dx « ApAx = 

(hN 1 / 3 /R 2 - h^^/R^iRx/N 1 / 3 - R2/N 1 / 3 ). The con- 
dition, Ri » R%,— > S/h = R1/R2. Therefore, the 
tunneling probability for a particle, T S i ng i e ~ e~ s '^ ~ 
e -fii/i? 2 \yhen all the particle, simultaneously, undergo 
such transitions, the compact star will tunnel through 
the energy barrier. For a star with particle number, 
N, the tunneling probability is, T w e -N(Ri/R 2 )^ p or 
typical white dwarf and neutron stars, particle number, 
N ~ 10 57 , which is a huge number, and therefore, the 
tunneling probability to the unstable ground state is ex- 
tremely small, and therefore, negligible. 



V. CONCLUSIONS 

We have analysed in details the effects of short range 
corrections to gravity on the compact stellar objects such 
as the white dwarf and neutron stars. The corrections 
lead to two new effects. One of them is the emergence of 
surface tension which is the characteristic of short range 
interactions and has been discussed in detail in our previ- 
ous publication [j| • The second effect is that the compact 
stars become metastable. We have shown that the tun- 
neling probability to the unstable ground state is so small 
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that the stars once trapped in the metastable state would 
remain there. 
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